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Abstract 

The excitation and absorption properties of grating coupled graphene surface plasmons were studied. It was found 
that whether a mode can be excited is mainly determined by the frequency of incident light and the duty ratio of 
gratings. In the structure consisting graphene bilayer, a blueshift of the excitation frequency existed when the 
distance between neighbor graphene layer were decreased gradually. In graphene-grating multilayer structures, a 
strong absorption (approximately 90% at maximum) was found in near-THz range. 
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Background 

Recently, a lot of work has been done based on graphene 
due to its unique properties in electric, magnetic, ther- 
mal, etc. [1-3]. Graphene is carbon atoms arranged in a 
two-dimensional honeycomb lattice, in which the elec- 
trons behave like massless Dirac fermions with linear 
dispersion [4,5]. Graphene has strong plasmonic effects 
which can be modified by gating, by doping, and so on 
[2]. A controllable optical absorption was also found in 
structured graphene [6,7]. 

Up to date, the graphene is modeled usually to be an 
extremely thin film with a conductivity cr, which consists 
of both intraband and interband from Kubo formula 
[7-9]. The intraband conductivity with Drude type plays 
a leading role when ho)//u c was small [10]. Both trans- 
verse electric (TE) and transverse magnetic (TM) have 
the dispersion relations at monolayer graphene with 
dielectric materials on two sides [10-12]. In other words, 
the charge carriers coupling to electromagnetic waves 
will produce a new surface wave, namely graphene 
surface plasmons (GSPs). 

In the previous works, many numerical approaches 
were used to study the structured graphene, for example 
the finite element method (FEM) [13], finite difference 
time domain (FDTD) [14], and others [6,15]. A strong 
plasmonic response of graphene has been demonstrated 
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in a square-wave grating with a flat graphene on top 
[15]. In which, the graphene-based plasmon response 
lead to a 45% optical absorption. In a periodic array of 
graphene ribbons, remarkably large GSPs result in 
prominent optical absorption peaks [13]. In multilayer 
graphene, the absorption spectrum can be decomposed 
into subcomponents [6], which is helpful in understand- 
ing the behavior of GSP coupling. 

In this paper, we studied the binary grating bounded 
by graphene on both sides. The rigorous coupled wave 
analysis (RCWA) [16,17] was used the first time as we 
know to characterize the graphene-containing periodic 
structures. The excitation condition and excitation in- 
tensity seemed to be influenced by the grating constant, 
duty ratio and the distance between the graphene layers. 
When introducing more graphene layers into the struc- 
ture periodically, a strong absorption band was found in 
the near-THz range. 

Methods 

Electromagnetic mode of binary grating-graphene 

Previous research has shown that the conductivity of 
graphene came from the contribution of intraband and 
interband [18-22]. The interband conductivity tends to 
be ignorable when hco < ft c (see [10]). Then the intra- 
band conductivity can be expressed as [23] 
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where \a c is the chemical potential, relating to the 
electron density. Equation 1 became a Drude type when 
fi c /k B T» 1, i.e. 



(2) 



First, we studied graphene layer positioned at the inter- 
face of two half planes with permittivity Si and e 2 , respect- 
ively. The electric induced current on graphene layer 
resulted in a magnetic field difference, which led to the 
coupled GSP on graphene layer. Using Maxwell equation 
and boundary condition, GSP modes were proved to 
existed for both TE and TM polarization [12,23-25]. For 
TE mode, the dispersion relation was as follows: 



/p 2 -k 2 0 e 2 + yj fi-klsi-icdfaG = 0, 
and for TM mode it became 
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Because the imaginary part of conductivity (2) was 
positive, no solution of Equation 3 was found in real, 
which meant the TE mode GSP could not be excited. 

For TM mode, put Equation 2 into Equation 4, we 
found 
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Here, we defined # eff = f}/k 0 = fic/co as the effective index 
of GSP. After making a transformation of (co, # e ff) —> ft), 





Figure 2 Binary grating graphene structures, (a) The bilayer 
graphene structure, (b) The multilayer graphene structure, h is the 
grating layer thickness. A is the period of grating. Li is the width of 
dielectric with L(L 2 ) is the width of dielectric with e 2 . The duty ratio 
is f 2 = L/A, and f y = 1 - f 2 . In this paper, we simply set ^ = 1 and e 2 = 4. 



the dispersion relations were obtained and plotted in 
Figure 1. The wave vector was normalized by k A0 = 2nl 
A 0 j A 0 = 1 |^m. As a local mode, GSP modes were same 
as the surface plasmon polaritons (SPPs). They cannot 
be excited directly from the air. And in our work, grat- 
ings were used to provide an external wave vector to 
match the phase condition. 




Figure 1 Dispersion relations of graphene surface plasmons (GSPs) 
on monolayer graphene with different material on two sides. 

Here, we use the graphene parameters of jj c = 0.2 eV, r" 1 = 1 meV. 
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Figure 3 Attenuated total reflection of the structure in Figure 

2a. A = 1 1 Mm, L = 1 0 urn, so the duty ratio is 10/11. Each of the 
absorption peaks (on blue solid line) corresponds to a GSP mode. 
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Table 1 The resonant frequency of different orders 



Order of GSP (N) 


1 


2 


3 


4 


5 


6 


7 


(jJq (meV) (RCWA) 


11.9 


16.7 


20.5 


23.7 


26.3 


28.9 


31.1 ... 


oj-i (meV) 


11.70 


16.61 


20.38 


23.55 


26.34 


28.86 


31.18 ... 


(theoretical) 

















w 0 was the numerical results obtained by RCWA. was the theoretical results 
from Equation 10. 



Rigorous coupled wave analysis in graphene-containing 
structures 

In Figure 2a, we used h to be the depth of grating (thick- 
ness of gratings). The h was also the distance between two 
graphene layers. In multilayer structures of Figure 2b, 2 h 
was the longitudinal period. The structures were designed 
to only contain two kinds of interfaces. 

In common, the conventional RCWA based on the 
Floquet's theorem [26] was unable to be used for the 
graphene-containing structures as the electric field will 
induce a current with current density J = <xE, while gra- 
phene was included. 

In RCWA, the field was expanded into the form of 

E x = Y, E M^ nX - at \ H y = ^Hyinyi^). 



(6) 



So the current density J can also be expended to the 
sum of spatial harmonics with different wave vector 
components. To obtain the reflection, transmission, ab- 
sorption, field distribution, and other optical properties 
of such structures as shown in Figure 2, a nonzero item 
must be included in the boundary condition of H y field 
considering the induced current, 



Hy+ —Hy 



(7) 



According to the principle of superposition, H y will 
also be continuous at the interface if each spatial har- 
monics subcomponent satisfied the boundary conditions 
independently, 
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Figure 5 Absorption spectra for duty ratio vs the frequency 
fixing the light path of grating period. 



Hy(n,x,y Q+ z)-H y (n,x,y Q _z) = oE z {n,x,y 0 z), (8) 

in which n was the order, ± in subscripts represented ap- 
proaching to y 0 from two different directions. After the 
modification on the RCWA program, we can utilize it to 
deal with the graphene-containing structures. 

Results and discussion 

Phase matching condition 

For a structure with a binary grating bounded by gra- 
phene layers on two sides shown in Figure 2a, the atten- 
uated total reflection spectrum is plotted in Figure 3 
using the modified RCWA method when it was illumi- 
nated normally. A set of absorption peaks each corre- 
sponding to a GSP mode was shown in blue solid line. 
From left to right, each peak corresponding to a GSP 
mode ordered with 1, 2 ... 

In the structure shown in Figure 2a, there exist two 
kinds of interfaces, i.e., ^-graphene-g! and ^-graphene-^. 
When GSP is propagating along the interfaces, the phase 
shifts were epi and cp 2 for the two kinds of interfaces, 




^ O loo ©e coco 

5pm I 



Figure 4 Field distributions of \E y \ of first (11.9 meV, left) and second (16.7 emV, middle) order GSP modes. The last figure was real part 
of E y of second order. 
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Figure 6 Corners of grating will become the scatting points of the incident light which was the source of GSPs. These scatting points 
can be divided into two kinds due to the geometric symmetry, which is A and B. Each scatting point will scatter into two GSP modes 
propagating in two directions (blue and green). 



respectively. 8 was the total phase loss considering two 
abrupt phase changes when GSP propagates across the 
joints between the two kinds of interfaces in a grating 
period. At the excitation frequency, the phase change in a 
grating period should satisfy the relation 



(9) 



which was known to be the phase matching conditions 
[27,28]. In Equation 9, N is the integer and can be 
rewritten as 



N(o)) 



1 



A 



(10) 



where f 2 = LI A and/! = 1 -f 2 , fi\ and f> 2 were the wave vec- 
tors of GSP on two kinds of interfaces, respectively. When 
N was a nonnegative integer, the GSP mode could be ex- 
cited, and N can be defined as the order of surface modes. 
The resonant frequencies can be obtained both from ab- 
sorption spectrum in Figure 3 and theoretically from 
Equation 10 (S = 0). They were given in Table 1 and agreed 
well for high order modes. But for low order modes, some 
deviations existed between numerical and theoretical 
caused by the coupling of GSPs on two graphene layers. 

The field distributions of orders 1 and 2 of the struc- 
ture in Figure 2a were given in Figure 4. It was indicated 
that the GSP field distributions had nodes as standing 
wave because the GSP modes propagating in two direc- 
tions were excited simultaneously. 

Duty ratio and stand wave interference 

By using the modified RCWA, the absorption spectrum 
was obtained in Figure 5 when varying f, where f=cp\I 
(epi + cp 2 ), <Pi, cp 2 had the same meaning as Equation 9. 
From the discussion above, when the phase match con- 
ditions were satisfied, GSPs could be coupled and ab- 
sorption peaks should appear. But some break points 
appeared as shown in Figure 5, which meant that GSPs 
cannot be excited at this frequency and this duty ratio 
and it was reasonable to believe that there must be other 
process related to /of controlling the excitation of GSP. 

From the field distributions in Figure 4, each corner of 
the grating was a singular point of field and these scatting 
points became the sources of surface wave, as Figure 6 
shown. In periodic, we only need to consider the scatting 



in one period, i.e., A and B. Each scatting point will couple 
to two GSP modes propagating in two directions. So, the 
field can be expressed in four terms, which is [28,29] 

p — p Q tf(kox-ci>i) _|_ p Qe i{-kox-at) _|_ p^ e i[k 0 (x-x 0 )-o)t+TT] 
_|_ p Qe i[-ko(x-x 0 )-6)t+7i] 



4F 0 sin KqX sin 



2 J 2 



(11) 



First two terms were GSP excited by one set of points 
(A in Figure 6) with two propagating directions (blue 
and green) and the last two terms were that from 
another set of points (B in Figure 6), where Xq IS the 
distance of A and B in the form of light path (koX 0 = 
= <Pi = (<Pi + ^2)/= 2nNf), Because in real space, dif- 
ferent interfaces (£i/£i and £1/82) had different propagat- 
ing constants, the expression might be complex. Here, 
the light path of x was used. It is found that scatting 
points A and B had a phase difference of 7T. This was 
caused by the different geometric symmetries. From 
Equation 11, when sin(/<o# 0 /2) = 0, i.e.,/= m/N ( m = 0, 1, 
A/), field amplitude F would always be 0, which meant 
that the field cannot be excited. It was a cancelation 
process of two sets of standing waves that are coherent. 
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Figure 8 The absorption spectrums of grating-graphene periodic 
multilayer structure. 'Layers', number of graphene layers, which is the 
odd number between 2 and 26. The frequency ranges from 0 to 
60 meV (approximately 14.5 THz). The figure inset is the reflections. 



So, for GSP mode of N, N + 1 of none absorption points 
appeared. 

Coupling of GSPs on different graphene layers and 
resonant frequency shift 

From Table 1, we can see that for higher order modes, the 
consistency between the theory and the numerical results 
from RCWA was better than that of the lower order 
modes. It was because the structure consists of bilayer of 
graphene and there could be interaction between GSP 
modes on neighbor graphene layers determined by the 
depth of the grating. In order to understand the behavior 
of GSPs coupling, in Figure 7, the absorption spectra 
were given as a function of the grating deepness /z. A 
blueshift of absorption peaks was found when the grat- 
ing became thin. The oscillator model is used to de- 
scribe this phenomenon of spectrum blueshift [30,31]. 



^ upper H~ Y% upper — Fq6 



-i{o)t-(p n ) 



-^0*upper-K(X A0)x\ o 



In Equation 12, k(w, /z, A0) was the coupling coeffi- 
cient and n, /z, and AO were order of GSP mode, thick- 
ness of grating, and phase difference of GSPs on two 
graphene layers, respectively. Essentially, the GSPs were 
surface waves so they interact with each other via evan- 
escent interactions, and the coupling intensity decayed 
exponentially with h increasing. For fixed /z, the lower 
order modes had larger skin depth (stronger coupling in- 
tensity) than the higher orders; then, the stronger coupling 
resulted in a large spectra shift. The phase difference of 
AO also had affection to the absorption frequencies. How- 
ever, in our case, the wavelength (15 meV~ 82.8 um) was 
much larger than the thickness of grating layer (h = 
10 um), it is reasonable to assume AO is approximately 0. 
This can also be obtained clearly from the field distribu- 
tion in Figure 4 that the electric fields on upper and 
lower graphene layers oscillated synchronously. This 
conclusion can still hold in multilayer graphene- 
grating structures. Finally, k{yi, /z, AO) * e~ hq(<n \ where 

q{n)=[p\n)-e 1 , 2 k^ 12 . 

Suppose the solution of having the form of x up = x<± own = 
x 0 e~ lcot (no phase difference between GSP on neighbor 
layers), it is found that the resonant frequency became 



Co 2 = o) 2 + K(h,n, AO)- 



(13) 



(12) 



When h was small (h < 4 um), the larger i<(n, h, AO) * e~ h 
was the larger shift of resonant frequency would be. 
And obviously, i<{n, h, AO) was approaching 0 rapidly 
when h was large enough, which meant that the reson- 
ant frequency became a stable value of (col-y 2 /l) 1 ^ 2 . 
Otherwise, i<{n, /z, A0) was also related to the order of 
GSP. The high order mode had a small skin deep with 
weak coupling intensity and less blueshift. When h 
tends to be 0, the grating became too thin to excite the 
surface mode. This was why the absorption disappeared 
when h = 0 in Figure 7. 




2.5Mm 



Figure 9 Field distributions. The real part (a) and (b) and magnitude (c) of E y in multilayer structure of different orders, (a) Excitation at the 
frequency of 24.6 meV. (b) and (c) Excitation at the frequency of 28.4 meV. 
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Strong absorption in grating-graphene multilayers 

Moreover, the behavior of multilayer structures shown 
in Figure 2b was also investigated using the modified 
RCWA and the absorption and reflection spectra were 
given in Figure 8. When increasing the number of gra- 
phene layers, it can be seen that the resonant frequen- 
cies do not change but for several lower order modes. 
Though the reflections were always weak within the res- 
onant range, it is obvious that the more graphene layers 
included, the stronger the absorption is (almost 90% 
when it contained 26 graphene layers). 

The field distributions of Figure 9 also give the same 
conclusion that the stand waves on each graphene layer 
were almost oscillated synchronously. The energy was 
mainly located and absorbed by the graphene layer as 
we expected. 

Conclusions 

On conclusion, the rigorous coupled wave analysis was 
modified to compute the excitation of graphene surface 
plasmon in graphene -containing binary gratings struc- 
tures. Under the phase matching conditions, the excita- 
tion of the graphene surface plasmonics was determined 
by the distance between graphene layers and duty ratio 
of gratings, and the mode suppression can be realized by 
modifying the grating constant and duty ratio. A blue- 
shift of the excitation frequency was obtained for en- 
hanced coupling between GSP of neighbor graphene 
layers. Increasing the number of graphene layers had 
almost no effect on the excitation frequency of GSP but 
would lead to a high absorption with negligible reflec- 
tion in near-THz range. Finally, the resonant frequency 
and absorptions can be easily modified by manipulating 
the structure parameter, including grating constant, duty 
ratio, and distance between the graphene layers and 
number of grating, and graphene-containing grating 
might become potential applications of THz region, such 
as optical absorption devices, optical nonlinear, optical 
enhancement, and so on. 
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